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Abstract 

Numerical approximate computation can solve large and complex 
problems fast. It has the advantage of high efficiency. However it 
only gives approximate results, whereas we need exact results in many 
fields. There is a gap between approximate computation and exact 
results. In this paper, we build a bridge by which exact results can be 
obtained by numerical approximate computations. 

Key words: Numerical approximate computation, Symbolic-Numerical 
computation. Continued fraction. 



1 Introduction 



Historically, the fields of symbolic computation and numerical computation 
have been developed by two distinct groups of people, having relatively lit- 
tle interaction and overlaps with each other. Symbolic computations are 
principally exact and stable. However they have a high complexity. There- 
fore, they are slow and in practice, are applicable only to small systems. 
Numerical approximate computation has the advantage of being fast, flex- 
ible in accuracy and being applicable to large scale problems. In Recent 
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two decades, numerical methods are applied in the field of symbolic com- 
putations. In 1985, Kaltofen presented an algorithm for performing the 
absolute irreducible factorization, and suggested to perform his algorithm 
by floating-point numbers, then the factor obtained is an approximate one. 
After then, numerical methods have been studied to get approximate factors 
of a polynomial [3] [10] |14j |15j |16j|17j . In the meantime, numerical methods 
are applied to get approximate greatest common divisors of approximate 
polynomials [I] [2] [II] [5j , to compute functional decompositions [6j, to test 
primality[9] and to find zeroes of a polynomial [12]. In 2000, Corless et al. 
applied numerical method in implicitization of parametric curves, surfaces 
and hypersurfaces[3]. The resulting implicit equation is still an approximate 
one. 

There is a gap between approximate computations and exact results [T8]. 
People usually use rational number computations to override the gap [7]. In 
fact, these are not approximate computations but big number computations, 
which are also exact computations. In 2005, Zhang et al proposed an al- 
gorithm to get exact factors of a multivariate polynomial by approximate 
computation [19] but they did not discuss how to override the gap. Com- 
mand convert in maple can obtain an approximate rational number from 
a float if we set variable Digits to a positive integer. However, in order to 
obtain exact rational number from its approximation, we need to know two 
things. One is at what accuracy the float should be obtained by numerical 
method; another is when we should stop and return the rational number 
we want by the continued fraction method. So we can not obtain the exact 
rational number from its approximation by command convert and variable 
Digits. In this paper, we solve the two problems, which can be described as 
follows: 

There is an unknown rational number m/n we want to obtain, and as- 
sume that there are approximate methods to obtain its approximation at 
arbitrary accuracy. We also know an upper bound of absolute value of 
its denominator in advance. The two problems will be solved such as: At 
first, we discuss how to determine e which is a function in N. And then use 
the approximate methods to obtain a floating-point number x, an approx- 
imation of m/n at accuracy e > 0, i.e. \x — m/n\ < e. Second, we give a 
criteria to stop our program and return the exact rational number we want 
by continued fraction method. 

The remainder of the paper is organized as follows. Section 2 gives a 
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review of continued fraction. Section 3 discusses how small the error needs 
to ensure the exact number to be obtained and how to get the exact number 
from its approximation. Section 4 gives some experimental results. The 
final section makes conclusions. 



2 Continued fraction 



A continued fraction representation of a real number is one of the forms: 

1 

aoH ■ 1 , (1) 



where oq is an integer and ai, 02, as, ■ • • are positive integers. One can ab- 
breviate the above continued fraction as [ao; ai,a2, ■ ■ ■]■ For finite continued 
fractions, note that 

[ao;ai,a2, • • • ,0^, 1] = [ao;ai,a2, ■ • • ,0^ + 1]. 

So, for every finite continued fraction, there is another finite continued frac- 
tion that represents the same number. Every finite continued fraction is 
rational number and every rational number can be represented in precisely 
two different ways as a finite continued fraction. The other representation is 
one element shorter, and the final term must be greater than 1 unless there 
is only one element. However, every infinite continued fraction is irrational, 
and every irrational number can be represented in precisely one way as an 
infinite continued fraction. An infinite continued fraction representation for 
an irrational number is mainly useful because its initial segments provide 
excellent rational approximations to the number. These rational numbers 
are called the convergents of the continued fraction. Even-numbered con- 
vergents are smaller than the original number, while odd-numbered ones are 
bigger. If successive convergents are found, with numerators /ii, /12, ■ • •, and 
denominators ki,k2, - ■ then the relevant recursive relation is: 

The successive convergents are given by the formula 

hn _ O-nhn-l + ^n~2 
kn dnkji—l + kji—2 

where = 1, /i_2 = 0, = and k-2 = 1. Here are some useful 
theorems^: 
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Theorem 1 For any positive x R, it holds that 

xkji—i + kfi^i 



[ao,ai, • • • ,an-i,x\ 



(2) 



Theorem 2 The convergents of [oq, ai, 02, • • •] are given by 

r 1 

[ao,ai, - ■ ■ ,an\ — -r- 

and 

knhn—l ~ kn—lhn = • 



Theorem 3 



and 



kri 



k 



1 



n-1 



< 



k'ri k 



n'^n—l 



< 



+1 



In order to recover exact rational number, we introduce a controlling error 
into the conventional continued fraction method. The continued fraction 
method is modified as follows. 



Algorithm 1 Continued fraction method 

Input: a nonnegative floating-point number a and e > 0; 

Output: a rational number b. 



Step 1: i := 1 and xi := a; 

Step 2: Getting integral part of Xi and assigning it to Ui, assigning 
its remains to bi. If bi < e, then goto Step\^ 

Step 3: i :=i + 1; 

Step 4: Xi := and goto Step\^ 

Step 5: Computing expression ([7P and assigning it to b. 

Step 6: return b. 



We will discuss the controlling error e in algorithm [T] in the next section. 
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3 Recovering the exact number from its approxi- 
mation 



In this section, we will solve such a problem: for a given floating number w 
which is an approximation of rational number ^ , how do we obtain integer 
m and n? Without loss of generality, we always assume that m,n are positive 
number. At first, we have a lemma as follows: 

Lemma 1 m, n,p, q are integer, and pn >0. If\^~p\<^! then ^ = p- 

Proof. 1^ — ^1 = ^^"^""^"^ . Noticing \pm — qn\ is a nonnegative integer, and 
\^ — p\ < ^5 yields \pm — qn\ < 1. Hence \pm — qn\ = 0. That is ^ = ^. 
The proof is finished. 

Corollary 1 m, n,p, q are integers andp > 0,n > 0. Let N > max{p, n, 2}. 
- p \ < N{N-1)' ^ ^ p- 

Proof: When p / n, it holds that pn < N{N - 1). Hence |^ - ^| < 
< . According to lemma [H it is obtained that ^ = f • When 
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N{N-1) - pn- ^^^wJ-^^i^B ic±ii±iia,LU lu lo wuLa,iiit>u uiiau „ — p- 

p = n, we have 

m q \m — q\ 1 n 1 

p' 77 N{N-1) ' ^' N{N-1)-N-1- 

So, m = q. The proof of the corollary is finished. 

And now, we study how small the error l^i' — -^1 needs so as to get ex- 
act rational number ^ from its approximation w . The following theorem 
answers this question. 

Theorem 4 Let x = ^ be a reduced proper fraction, and N > max{n,2}. 
Assume that \x — w\ < \/{2N{N — 1)). // we get positive rational number 
p/q such that \p/q — w\ < 1/{2N(N — 1)) , where q < N , then it holds that 
X = q/p. 

Proof: From the assumption of the theorem, we have \x — q/p\ < 1/[N{N — 
1)). According to corollary [H it holds that q/p = m/n = x. The proof of 
the theorem is finished. 
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We have answered how smah the error is so that we can recover the exact 
rational number from its approximation. The remaining problem is how to 
get the exact number. We attack it by continued fraction method. 

Let n2/ni be a rational number and tq its approximation. Their continued 
fraction representations are n^jnx = [ao, ai, • • • , ol] and ro = [6o> ^i, • ■ ■ ; 
respectively. We wish that Oj = ftj for i = 0, 1, 2, ■ • • , L — 1 and for the last 
term of the continued fraction representations of n^jnx^ either a/, = hi, or 
a/, — 1 = 6i, so that we can get n2/ni from [6o, ^i, • ' ' ; ^L+i]- This is the 
following theorem: 



Theorem 5 Let nijux he a rational numher and its approximation. As- 
sume that n2,ni are coprime positive numbers, where n2 < ni,and ni > 1. 
The representations of n2/ni and r^ are [oq, ai, • " " ) cll] o-nd [hQ,hi, ■ ■ ■ , 6^./] 
respectively. If |ro — n2/ni\ < l/(4ni(ni — 1)), then one of the following 
statements must hold. 

• ai = bi (i = 0,l,--- ,L); 

• ai = bi (i = 0,1, - ■ ■ ,L - 1), ai - 1 = hL, and hL+i = 1. 



According to assumption of n2 < we have that oq = 0, and bo = 0. Hence 
oq = bo- In order to finish the proof of theorem [H we need two lemmas. Due 
to n2/ni = [ao, ai, • • • , ql] and tq = [ho, 5i, • • • , 5^], we have the following 
expansions: 

ni ns 712 n4 ul-i , 1 

— = Ol H , — = 02 H , • ■ ■ , = ai-i H , riL = ol (3) 

n2 n2 ns hl ni 

and 

— = bi + ri, — = 62 + r2, • • • , =hL+rL, ■■■ , = hM (4) 

Denoting di = ri — n^+2/7^^+l, we have a lemma as follows: 



Lemma 2 Let 712/ ni he a rational numher and ro its approximation. As- 
sume that 712, ni are coprime positive integers, where n2 < ni,and ni > 1. 
The representations of n2/ni and ro are [ao, ai, ■ ■ ■ , ai] and [60, ^i, ■ " " ; bu] 
respectively. And assume that Oj = hi for i < k < L(k is a positive integer). 

Then when [dk[ < — ^tt; holds that a^+i = for k < L — 1; 

when < „ 1^ I -|\ , it holds that ai = h^ or ql — I = b^. 
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Proof: At first, we show that under the assumption of the lemma if we have 

1 



nk+2{nk+2 + nk+idk) 



< 



nk+2 



(5) 



then, it holds that a^+i = b^+i for < L — 1, and ak+i = or a^+i — 1 = 
bk+i for k = L — 1. We discuss it in two cases: 

Case l{k < L—1): From = rk—nk+2/nk+i, it holds that = nk+2/nk+i+ 
d^. Hence we have that 



1 



nk+1 
nk+2 



n 



k+l 



dk 



1 



nk+i 



n 



k+l 



dk 



nk+2{nk+2 + rik+idk) nk+2 nk+2{nk+2 + nk+idk) 



1 

rk 



o-k+i + 



nk+3 



n 



k+l 



dk 



nk+2 nk+2{nk+2 + nk+idk) 
Hence, it is obvious that Uk+i = bk+i if and only if 

nk+3 n\ 



bk+i + Tk+i 



< 



^k+l^k 



nk+2 nk+2{nk+2 + nk+idk) 



< 1. 



(6) 



Therefore, if inequality ^ holds, then above inequality is guaranteed. 
Case 2: (when k = L — 1) We have 



1 



aL 



nldL- 



1 



ri-i nL+i{nL+i + nidi-i] 

From the above equation, if 



bl + TL 



nidL- 



ni+iini+i + nidi-i] 



< 



1 



nL+2 



, then aL = b^ for di^i < 0, and — 1 = b^ for d^^i > 0. Therefore, we 
have shown that if inequality ([5]) holds, then Uk+i = bk+i for k < L—1, and 
afc+i = bk+i or ak+i - 1 = bk+i for A; = L - 1. 



On the other hand, we have 



n 



k+l 



dk 



n 



k+i\dk\ 



nk+2{nk+2 + nk+idk) 
nk+i\dk\ 



< 



\nk+2{nk+2 + nk+idk)\ 
nk+i\dk\ 



nfc+2(|;^ + 4|) nk+2\ 



' nk+2 



141)1 
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So, in order to ensure inequality ([5]), we only need it holds that 

nk+i\dk\ ^ 1 



nfc+2|(^ - 141)1 nk+2 
Solving inequality ([7]) yields 



I , I ^ nk+2 

nk+i[nk+i + 1) 



When < L — 1, we have that > 1- So, it holds that 

1 ^ 



nfc+i(nfc+i - 1) nfc+i(nfc+i + 1) 
Accordingly, it is obtained that 

141 < (9) 

nfc+i("-fe+i - Ij 

When k = L — 1, we have that n^+i = 1, so it is obtained that 

|4-ll<^^— -T (10) 

nL[nL + 1) 

The proof of lemma [2] is finished. 



Lemma 3 Letn2/ni be a rational number and r^ its approximation, where 
712, ni are coprime positive integers, andn2 < ni,andni > 1. The continued 
fraction representations ofn2/ni and x are [ao, ai, • • • , a^] and \bQ,bi, ■ ■ ■ ,bM] 
respectively. Denote di = r j — nj-(_2/'^i+i for i = 0, ■ ■ ■ , L. Assume that ai = 
h for i<k<L—l(kisa positive integer ). Then when [41 < nk+\{nk+i-i) ' 
it holds that 

|4+il < 7 7t|4| (11) 

nk+2{nk+2 - 1) 

Proof: Under the assumption that Oj = bi for i = 0,1,-- - ,k, from equation 

(El), we get 4+1 = / — j-T- Hence we deduce a relation as 

follows: 

,2 



%+il4l _ nfc+i|4 



'■k+2 



+ nk+ink+2dk nfc+2(;^ + 4) 
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nk+i{nk+i - l)\dk\ 



nk+2 r^t':'-'^ + K+i - 1)4) 



nk+i{nk+i - l)\dk\ 



nk+2{nk+2 - 1 + ^^Hgr^ + K+i - ^)dk) 



When 141 < ^, it holds that "fc+i "w ,^/^ i)^ > q. Hence 

we have a relation between 4+i and 4: 

14+1 [ < 7 tt|4| 

nk+2[nk+2 - IJ 

The proof of the lemma is finished. 

And now, let us prove the theorem. If |4| = \fo — n2/ni\ < l/(4ni(ni — 
1)), From lemmaO we can get \di\ < l/(4nj_|_i(nj_|_i — 1)) for i = 0, • • • , L — 1. 
Note that ni > ni+i = 1 and 

1 1 1 

< 7 7T < 



4ni+i(ni+i - 1) nj+i(ni+i + l) nj+i(nj+i - 1) 
when rii+i > 1. So, it holds that 

\di\ < -, — < 



4ni+i(ni+i - 1) ni+i(ni+i + l) 

for i = 0, ■ ■ ■ , L — 1. According to lemma [21 the proof of the theorem is 
finished. 

For an unknown rational number n2/ni and its approximation ro, theorem 
[5] shows that n2/ni = [bo, • • • , 6l] or n2/ni = [bo, bi, - ■ ■ , bL, 1] when |ro — 
n2/ni[ < l/(4ni(ni — 1))- However, we do not know what the number L 
is. If we make large enough when n2/ni = [6o, • ' ' j ^1/]; or make 6l+2 
large enough when n2/ni = [6o, 6i, ■ • • , bi, 1], then we recover n2/ni easily. 
The following theorem solve this problem. 

Theorem 6 Let 77-2/^1 be a rational number and ro its approximation. As- 
sume that n2,ni are coprime positive integers, where n2 < ni,and ui > 1. K 
is a positive integer. The continued fraction representations ofn2/ni andr^ 
are [oq, ai, • • • , a/,] and \j3o,bi, ■ ■ ■ ,bM] respectively. //|4| = — "'2/?^i| < 
1/{{2K + 2)ni(ni — 1)), then one of the following two statements must hold 

• Qi = bi for i = 0, ■ ■ ■ , L, and > K; 
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ai = bi for i = 0, - ■ ■ ,L - 1, and = aL - 1, bL+i = 1, 6l+2 > K. 



Proof: From equation ([3]) and equation dH), we have that ^ = a^-i + 
riL = ai and = + When |do| = \ro-n2/ni\ < 1/{{2K + 

2)ni(ni-l)) < l/(4ni(ni-l)), fromlemmaEl it holds that \di\ < 1/{{2K+ 
2)ni+i(ni+i — 1)) for i = 0, 1, • • • , L — 1. Furthermore, from theorem [5] it 
holds that ul-i = bi-i and di-i = vl-i — Tr~- So we have 

1 _ _ nldL-1 

TL-l 1 + nLdL-1 



We discuss it in two cases: 



Case 1 (di-i > 0): Due to — j^^^^j^ < 0, we have 
1 - . ... nldL-1 . 



TL-l 



aL - 1 + (1 



1 + ULdL^ 



bL + rL 



Hence 



rL = (l 



nldL-1 X _ 1 + nLdL-i - n|dL_ 



1 + ULdL-i 1 + nLdL-i 

^ J_ _ 1 + nLdL-i 

rL 1 + nLdL-i - nldL-1 

It is obvious that 1 < l/r^ < 2 when [d/^-i] < l/(4n/,(n/, — 1)). Therefore, 
we have 



1 + 



1 + riLdL- 



rL 1 + nLdL-i - nldL-1 

We get that = 1 and 

1 + TiLdL-i - nldL- 



1 = 1 + 



nldL-1 



1 + ULdL-i - nldL-1 



nldL-1 



bL+2 + ?^L+2- 



(12) 



So, 



1 + riLdL-i - nldL-1 
nldL-1 



where [.] stands for getting the integral part of a number. 
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We hope that 6^+2 is greater than some integer K, which is used as a 
sign that ho, - ■ ■ ^hiA have been obtained. Solve the following inequality: 

1 + nidL-i - n\dL-i ^j^^., , 2, ^7^2, 

— > K 4^ 1 + ULdL-i - n^dL^i > Knj^dL-i 



nidi 



^ di-i < 



1 



{K + l)nl-nL 
Since {2K + 2)nL{nL - I) > {K + l)n\ - ul, we take 



dL-i < 



1 



{2K + 2)nL{nL-l) 



Case 2{dL-i < 0): Due to -j^^j^ > 0, we have 



So, ai = bL and 



It is obtained that 



Therefore 



ai + (- 
rL = ( 



nldL-i 



1 

rL 



1 + riLdi^i 

n\dL-i 
1 + nidL-i 

1 + riLdL- 
n\dL-i 



hi + rL 



- = hi+i + r-L+i 



1 + riLdi- 
n\dL-i 



Solve the following inequality 

1 + TlLdl-l 



nldi-i 
^ i-di-i) < 



> K ^1 + nidL-i > nU-dL-i)K 



1 



Knl + riL 



It is obvious that {2K + 2)nL{ni — 1) > Kn\ + for ni > 2. So we take 



-dL-i < 



1 



{2K + 2)nL{nL-l) 
The proof of theorem [6] is finished. 

For practical purpose, we hope the restriction on ni > 1 and ni > 712 can 
be lifted. So we have following theorem: 
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Theorem 7 Let no/^i he a reduced rational number and r its approxima- 
tion. Assume that nQ,ni are positive integers and N > max{ni,2}. K is 
a positive integer. The continued fraction representations of no/ni and r 
are [oq, ai, • • • , a/,] and [60, • • • , &Af ] respectively. If \d\ = |r — no/ni| < 
1/{{2K + 2)N{N — 1)), then one of the following two statements must hold 

• Gi = bi for i = 0, ■ ■ ■ , L, and > K; 

• ai = hi for i = 0, - ■ ■ ,L - 1, and bL = ai- 1, = 1, 6l+2 > K- 

Proof: We prove the theorem in three cases: 
Case 1 (ni > 1, no < ni): This is theorem [6l 

Case 2(ni = 1): We have that ao = uq/ui. li d = r — uq/tii > 0, then 
60 = ao andro < 1/((2E: + 2)2(2- 1)). So, it holds that 1/ro > (2(2A" + 2)). 
Therefore, it is obtained that hi > K. If d = r — n^/ni < 0, then 

|d| ^ ro = 1 - \d\ 



1 



So, we have that 60 = ^0 — 1, ^1 = 1 and 62 > K. 

Case 3(no > ni): From no/^i = ao+n2/ni, it holds that no/ni— oq = n2/ni. 
On the other hand, we have that |no/ni-r| < 1/{{2K+2)N{N -I)) < l/ui. 
So, we can deduce that oq < t < + 1. Accordingly, it holds that 69 = ^o- 
Hence, we have 

\d\ = |r - no/ni| = |6o + ro - ao - n2/ni| = |ro - n2/ni| 
= do< 1/{{2K + 2)N{N - 1)) < l/((2i^ + 2)ni(ni - 1)) 

Since ni > 1 and 712 < ni, from theorem [UJ the theorem holds. Therefore, 
the proof is finished. 

From theorem [71 we can get exact non- negative number 77-2 / ni from its 
approximation r as follows: 

Step 1: estimating an upper bound of the denominator of n2/ni. De- 
noted by N] 



r = ao - |d| = ao - 1 + 1 - = 60 + 1 

\d\ . \d\ 



1/ro = 1 + 



bi + 



l-\d\ 



i/n 



1 



\d\ \d\ 
62 > i/n - 1 > 4K + 2 > K 



1 > (4K + 4) 
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Step 2: computing 

^ ~ (2A^ + 2)iV(iV - 1)) 

Step 3: obtaining r by approximate method such that \r — n2/ni\ < d; 

Step 4: taking e = 1/N in algorithm [1] and calhng algorithm [1] to get 
b. So 6 = n2/ni. 

4 Experimental results 

The following examples run in the platform of Maple 10 and PIV 3.0G, 512M 
RAM. They take little time for obtaining exact rational numbers from their 
approximations, so we do not show time. 

Example 1. Let a be unknown rational number. We only know a bound 
of its denominator N = 170. According to theorem [TJ Computing rational 
number a as follows: Take K = 170+1, d = {2*K+2)*N*{N-1) = 9883120, 
and compute 1/d = 1/9883120. Assume that we use some numerical method 
to get an approximation b = 0.8106507864 such that |a — 6| < 1/d. Taking 
e = 1/K, we recover number a by algorithm[Tl We get [0, 1,4,3, 1, 1,4,] by 
the first 7 steps. When doing at step 8, we get 314, which is larger that K. 
We stop and return a = [0,1,4,3,1,1,4]. It is 137/169. 

Example 2. Let a be unknown rational number. We only know a bound 
of its denominator = 1790. According to theorem [71 Computing rational 
number a as follows: Take K = 1790 + 1, yields d={2*K + 2)*N*{N- 
1) = 11477079040. Assume that we use some numerical method to get an 
approximation b = 0.1788708777 such that |a — 6| < 1/d. Taking e = 1/K, 
we recover number a by algorithm [H We get [0,5,1,1,2,3,1,6,2,] by the 
first 8 steps. When doing at step 9, we get 2722, which is larger that K. 
We stop and return a = [0,5,1,1,2,3,1,6,2]. It is 320/1789. 

Example 3. Let a be unknown rational number. We only know a bound 
of its denominator N = 18. According to theorem [71 Computing rational 
number a as follows: Take K = N + 1, yields d = {2*K + 2)*N *{N -1) = 
12240. Assume that we use some numerical method to get an approximation 
1.882434634 such that [a — 6| < 1/d. Taking e = 1/K, we recover number a 
by algorithm [H We get [1, 1, 7, 1, 1] by the first 5 steps. When doing at step 
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6, we get 41, which is larger that K. We stop and return a = [1, 1, 7, 1, 1]. 
It is 32/17. 

Example 4. This example is an application in obtaining exact factors 
from their approximations. Let p = —16 — 56 * y — 48 * z + 64 * — 32 * 
X *y -\- * X * z — 45 — S)Q*y * z — 27 * he sl polynomial. We want to 
use approximate method to get its exact factors over rational number field. 
First, we transform p to a monic polynomial as follows: 

2 1 3 45 2 3 27 2 7 3 1 

p = x xy + -xz -y yz -z y — -z — - 

^ 2^ 4 64^ 2^ 64 8^ 4 4 

the least common multiple of denominators of coefficients of polynomial 
p{x,y,z) is 64, which is an upper bound[20] of denominators of coefficients 
of the monic factors of polynomial p. taking K = 66 yields d = (2 * K -\- 
2) * 65 * 64 = 557440 and e = 1/K. We use numerical methods to get its 
approximate factors as follows |20j: 

gi = 1.000000000000x+.6250000000067y+1.124999999530z+.5000000000000 

52 = 1.000000000000x-1.125000000015y-.3749999995480z-.5000000000000 

the error of coefficients of gi and g2 is less than 1 /d by the numerical meth- 
ods. According to theorem [TJ taking e in algorithm [U we obtain two exact 
factors: 

9 3 1 

02 = X y z 

^ 8 8 2 

5 Conclusion 

This paper builds a bridge spanning the gap between approximate computa- 
tion and exact results. The exact results can be obtained by our algorithm 
as long as we get a bound N of absolute values of their denominators and 
their approximations with a error less than 1/((2N + 2)N{N — 1)). Basing 
on our algorithm, we have succeed in obtaining exact factors of polynomi- 
als from their approximate factors. Our method can be applied in many 
aspect, such as proving inequality statements and equality statements, and 
computing resultants, etc. Thus we can take fully advantage of approximate 
methods to solve larger scale symbolic computation problems. 
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